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Abstract 



Fission events from Special Nuclear Material (SNM), such as highly enriched uranium or plutonium, can produce 
simultaneous emission of multiple neutrons and high energy gamma-rays. The observation of time correlations between 
any of these particles is a significant indicator of the presence of fissionable material. Cosmogenic processes can also 

^1 mimic these types of correlated signals. However, if the background is sufficiently low and fully characterized, significant 
changes in the correlated event rate in the presence of a target of interest constitutes a robust signature of the presence of 
£^*SNM. Since fission emissions are isotropic, adequate sensitivity to these multiplicities requires a high efficiency detector 
with a large solid angle with respect to the target. Water Cerenkov detectors are a cost-effective choice when large solid 
angle coverage is required. In order to characterize the neutron detection performance of large-scale water Cerenkov 

t-H detectors, we have designed and built a 3.5 kL water Cerenkov-based gamma-ray and neutron detector, and modeled 
the detector response in Geant4 pQ. We report the position-dependent neutron detection efficiency and energy response 
of the detector, as well as the basic characteristics of the simulation. 
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1. Introduction 

Legitimate cross border trade involves the transport 
of an enormous number of cargo containers. In order 
to verify that these containers are not transporting SNM 
without impeding legitimate trade, there is a need for 
fast, highly efficient, and large detectors that are rela- 
tively inexpensive. Such detectors must produce consis- 
tent yet distinct responses to both SNM and background, 
so that their effectiveness is not reduced by false positive 
or negative detections. They also need to have limited 
sensitivity to background radiation, such as cosmic ray 
induced background or Naturally Occurring Radioactive 
Material (NORM) present in certain legitimate forms of 
cargo. Both of these may contribute to false positives or 
reduce sensitivity to real SNM. 

SNM can either spontaneously fission or be induced to 
do so by an external source of gamma rays or neutrons. 
Since cargo containers are large, they can contain a sig- 
nificant amount of shielding. The fission emissions most 
likely to penetrate the container and interact with a de- 
tector are neutrons or high energy (greater than 3 MeV) 
gamma-rays. We propose that a water Cerenkov detector 
doped with a neutron capturing agent (such as GdCl3 salt) 
would be ideal for this application. Such a detector, sen- 
sitive to short timescale correlations between events, has 
a number of advantages; it is relatively inexpensive, non- 
flammable and noncombustible, environmentally safe, and 
easy to operate. 

Cerenkov detectors do not produce correlated signals 
from single fast neutrons in the same way organic scin- 
tillator does: in scintillator, fast neutrons are capable of 



producing a correlated signal via proton recoil followed by 
neutron capture. Water Cerenkov detectors only use the 
thermal neutron capture and prompt gamma-ray signals. 

Thermal neutron capture on natural Gadolinium has 
an extremely high cross section (49,000 barns). On cap- 
ture, a gamma-ray shower with energies adding to approx- 
imately 8 MeV is produced, and Cerenkov radiation pro- 
duced by the resulting Compton scatters is detectable by 
ordinary PMTs. SNO and Super-Kamiokande have shown 
that the Cerenkov process can generate enough photons 
in water to detect neutron captures or gamma-rays with 
an energy of approximately 3 MeV or greater, so long as 
the photocathode coverage is high (~40%) [2j [3]. Our 
group has since demonstrated the viability of this tech- 
nique above ground, operating 250 liter water based neu- 
tron detector with a photocathode coverage of 10%; reflec- 
tive detector walls made a lower photocathode coverage 
possible [5]. 

2. Detector Description 

In order to measure the neutron detection efficiency 
and maximize the performance of a large-scale detector, we 
have constructed a 3.5 kL water Cerenkov-based neutron 
detector. The detector consists of a cylindrical polyethy- 
lene tank with two PMT arrays arranged on the bottom 
and top, each with twenty Hamamatsu R7081 10 inch 
PMTs: although the current generation of R7081 PMTs 
are high quantum efficiency, the model used here does not 
have that improvement. The total photocathode coverage 
is approximately 19%. The tank is approximately 1.5 me- 
ter high, with a 2 meter diameter. The detector sits atop 
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Figure 1: The detector schematic showing the port and top hat for 
calibration arm insertion. The detector volume is in-between two 
PMT arrays supported by acrylic rods, and the detector sits inside 
a secondary containment vessel to protect against accidental leaks. 
Lead bricks shield against background gamma radiation from the 
floor. 
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Figure 3: A Hamamatsu R7081 PMT with factory water proof pot- 
ting (a), and the bottom 20 PMTs installed in the bottom detector 
array (b). 



in contact with the water are therefore constructed of ei- 
ther plastic or glass. The PMTs are mounted on a frame 
constructed out of acrylic and white polypropylene that 
sits inside the tank, shown in Figure [3^i. Each detector 
quadrant contains ten PMTs, five each on the top and bot- 
tom, with the top PMT array supported by acrylic rods. 
Because the frame structure and PMTs are buoyant, stain- 
less steel bars (sealed inside polypropylene bags) are used 
as ballast. Figure [3}}, shows the bottom four quadrants 
with the support rods in place. To increase light collec- 
tion, the detector wall in between the two PMT arrays is 
lined with UV reflecting Teflon. 

The detector is filled from the bottom using a water 
purification system capable of obtaining ultra pure deion- 
ized water (resistivity greater than 17 M^-cm). Municipal 
water is passed through three large de-ionizing (DI) resin 
bottles and is then circulated through a purification sys- 
tem purchased from South Coast Water Inc., consisting 
of a UV sterilizer, 5 micron and 0.22 micron filters, and 
an additional DI unit. Each unit in the system can be by- 
passed through a series of valves. After purification, the DI 
unit is by-passed and the water is doped with GdCl3. Any 
particulate matter and remaining biological contamination 
from the GdCl3 is eliminated by continued re-circulation 
through the filters and UV sterilizer. Finally, the doped 
water is sent to the detector inlet. 



Figure 2: The neutron calibration arm: the PMT reads out gamma 
radiation coming from the AmBe source. 



a layer of lead bricks to shield against background gamma 
radiation from the floor. Figure [I] shows a schematic of 
the detector: the white cylindrical opening in the center, 
called the top hat, allows for deployment of calibration 
sources at several locations inside the detector through a 
10 cm diameter portal in the lid. The calibration source 
sits at the end of a long arm (the calibration arm), shown 
in Figure [2] and discussed further below. 

It has been demonstrated that GdCl3 doped water is 
not compatible with stainless steel, contact with which re- 
sults in reduced water clarity [4] . All detector components 



2.1. Data Acquisition 

The DAQ and trigger electronics are mostly commer- 
cial VME and NIM modules. The only modules designed 
and fabricated in-house are eight-channel signal pick-off 
modules. The pick-off modules are needed to separate the 
fast PMT signal from the high voltage bias, which is car- 
ried on the same cable coming out of the PMT. 

From there, the PMT signals are passed through a 
Mini-Circuits 15542 BBLP-39+, 23 MHz low-pass filter 
in order to stretch out the signal, then amplified and dig- 
itized using CAEN V975 amplifiers and 200 MHz Struck 
SIS3320 digitizers. There are two triggers for the system: 
an internally generated trigger for physics events in the 
water and an external trigger used for calibration. The in- 
ternal trigger is formed from a four-fold coincidence among 
a group of 16 bottom PMT channels. The PMT signals 
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are fed to a CAEN V814 discriminator, set to trigger at 
-1 photo-electron. A CAEN V1495 FPGA registers the 
4-fold coincidences and issues the trigger to the waveform 
digitizers. 

3. Calibration and Detection Efficiency 

Calibration of the detector is needed to equalize the 
PMT gains and to set an approximate energy scale for 
physics events. We have three artificial sources: an LED 
for PMT gain, a 252 Cf fission source for neutrons, and 
a tagged neutron source to measure the neutron detec- 
tion performance on an event-by-event basis. A polyethy- 
lene calibration arm was constructed in collaboration with 
a group at Harvey Mudd College to deploy the neutron 
sources at multiple positions inside the detector. 

The tagged neutron source consists of an americium 
beryllium (AmBe) source and a Scionix type 51B51/2M- 
El-BGO crystal and PMT deployed together inside the 
calibration arm. AmBe sources emit a neutron in coin- 
cidence with a 4.4 MeV gamma-ray. Detection of this 
gamma-ray by the BGO detector forms a "tag", indicat- 
ing that a neutron has been emitted from the source. A 
rendering of the crystal, PMT, and AmBe source inside 
the polyethylene arm is shown in Figure [2] The arm was 
designed so that it can reach several positions inside the 
detector, making position dependent neutron detection ef- 
ficiency measurements possible. 

3.1. 252 Cf Calibration 

Three datasets were taken with a 252 Cf source posi- 
tioned at various locations outside the detector: seven 
inches from the detector wall, one meter from the detector 
wall, and two meters from the detector wall. The first po- 
sition is used to calibrate the response from Monte Carlo 
and to establish the quality of event-level cuts on the data. 
The one and two meter positions are used to determine 
how the event-level cuts perform with distant sources. 

Since high energy events, such as muons traversing the 
detector tend to saturate the response of the PMTs, we 
first screen our data to remove all events that contain at 
least one saturated PMT. Event level cuts are then applied 
to select neutrons in the detector. Each is determined by 
maximizing a Quality Factor, Q, defined as the significance 
after the cut is applied divided by the significance before 
the cut is applied. For a given parameter distribution di- 
vided into n independent bins, and for an analysis cut at 
the jth bin, the quality factor approaching from the left, 
or lower bound, is 
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Figure 4: In (a), the correlated (red) and uncorrected (black) 
summed response with the 252 Cf source seven inches from the detec- 
tor edge: the high-energy tail is due to cosmic-ray induced muons 
and gamma-rays. The charge quality factor for left (red) and right 
(black) cuts is shown in (b), where the maximum in the quality factor 
is used to establish the cut values. 



where On and Off represent the value of the parameter 
with and without the neutron source present, respectively. 
Q is defined to have a value equal to 1 when no cut is 
applied, i.e. when j equals n. Once a lower bound cut value 
has been determined, we then apply the same formula in 
reverse, or from the right, to determine the best jth bin 
for an upper bound cut: 
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The quality analysis for charge is shown in Figure [4] 
Because the fission source increases the number of uncorre- 
cted gamma-rays interacting in the detector, using the on 
source uncorrelated background as the no source present 
histogram for the quality analysis provides better rejec- 
tion of gamma energies. In Figure [4^i, the correlated on 
source data is shown in red and the uncorrelated on source 
is shown in black. Figure shows Q for the left (pink) 
and right (blue) cuts. The resulting cuts are applied to 
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Figure 5: In (a), the inter-event time for events passing the charge 
cut (red), and the fit to the uncorrelated background (black). The 
time quality factor for left (red) and right (black) cuts is shown in 
(b), where the maximum in the quality factor is used to establish the 
cut values. 



both the current event and the previous event to establish 
neutron-neutron pair events in the detector. 

Since the presence of a neutron source results in an in- 
crease in the uncorrelated trigger rate, a comparison of the 
inter-event time distribution with and without the source 
is not appropriate. Instead, the random trigger rate in the 
on source data is fit to an exponential, and a histogram 
filled randomly from the fit values serves as the no source 
present histogram. The results shown in Figure [5] with the 
on source inter-event time, as well as the quality factor. 

In order to pick out neutron-neutron events, we cut 
on both the current event and the previous event charge. 
Finally, we reject events in which the time difference be- 
tween the current event and the last muon is less than 46 
/is, which acts as a muon veto; the value is chosen as the 
same maximum inter-event time value allowed. A muon 
veto cut of 46 [is rejects 9.9% of the correlated background 
and increases the dead- time of the detector by 2.5%. Al- 
though there is modest improvement in rejecting muon- 
induced backgrounds, the majority of cosmic backgrounds 
appear to result from muons in the vicinity of the detec- 
tor but not traversing it. The final cut values for all four 
event-level cuts are shown in Table [TJ 



Current Charge 16 pe 72 pe 

Previous Charge 16 pe 72 pe 

Inter-event Time 4 [is 46 [is 

Muon Veto >46 [is N/A 

Table 1: Analysis cuts obtained by maximizing the signal significance 
between the background data run and 252 Cf data run with the source 
seven inches from the detector. 

Event Rate (Hz) 252 Cf Source No Source 



Raw 

Singles 

Doubles 



2390 

920 

55.0 



1770 

468 

13.7 



Table 2: Average event rates over 20 seconds for 252 Cf data run with 
the source seven inches from the detector, compared with no source 
present. The singles rate has only the energy cut and a muon veto 
applied, and the doubles rate has all cuts applied. 



Table [2] shows the event rates averaged over 20 seconds 
for the data run with the 252 Cf source seven inches from 
the detector edge compared to no source present. At this 
position, the solid angle coverage is 28% of 47r. The singles 
rate, or the rate of single neutron events in the detector, is 
the rate after both an energy cut on the current event and 
the muon veto cut has been applied. The doubles rate, or 
the rate of neutron-neutron pair events, is the rate after 
all cuts in Table [I] have been applied. 

3.2. Neutron Detection Efficiency 

The neutron detection efficiency is determined using 
the calibration arm described above. The AmBe gamma- 
ray spectrum from the BGO crystal in Figure [6] shows the 
primary 4.4 MeV gamma-ray peak. The smaller peaks at 
3.9 MeV and 3.4 MeV are also due to 4.4 MeV gamma- 
rays, where pair production and subsequent positron an- 
nihilation results in either one or two 511 keV gamma- 
rays escaping from the crystal. To maximize the signal 
to background ratio of our tag, we select events in the 
range 3 MeV to 5 MeV, where the background gamma- 
ray rate in the crystal is very small. Figure [7] shows the 
timing distribution of delayed detector events after our se- 
lection of AmBe tags. This inter-event time distribution 
can be well parameterized by a sum of two exponentials, 
representing a correlated and uncorrelated set of events. 
The correlated set represents neutron capture events due 
to the AmBe source in our detector. The mean capture 
time is 35 [is. This is consistent with the expected cap- 
ture time of thermal neutrons in water doped with 0.1% 
gadolinium [5j EJ [3 [8] . The number of neutron captures 
detected can be estimated by subtracting the exponential 
fit to the un-correlated events. The efficiency is taken to 
be the integral of the correlated events, the black curve 
in Figure [7| divided by the total number of neutron tags 
(after accounting for the crystal's background rate), or the 
red curve in Figure [7[ 
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Figure 6: The charge from the BGO crystal with (red) and without 
(black) an AmBe source present. The energy cuts around the 4.4 
MeV gammas and two escape peaks are shown in blue. 



Figure 8: The neutron detection efficiency verses the AmBe source 
position. For the position outside the detector (125 cm), the effi- 
ciency has been scaled by the solid angle coverage of the detector. 
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Figure 7: The time between tag and detector triggers for all events 
(red). In the 300-800 /zs range are the uncorrelated accidentals. The 
correlated events in the 0-100 /is range, with the uncorrelated back- 
ground subtracted (black), indicate a capture time of 35 /is. 



Three AmBe data sets were taken inside the detector 
to determine the efficiency and energy response at various 
radii. We assume that the detector is radially symmetric. 
One dataset was taken outside the detector, for which the 
efficiency is multiplied by the fraction of solid angle cal- 
culated from the source position. The position-dependent 
efficiency is shown in Figure [8] it ranges from 69.9% at 
the center of the detector to 31.3% outside the detector. 
The efficiency drop is expected near the detector edge, as 
neutrons can leave before being captured and many of the 
neutron capture gamma-rays escape the detector before 
interacting. 

3.3. Monte Carlo Response 

A full Monte Carlo was written in Geant4 to model 
our detector. Our objective is to utilize the Monte Carlo, 
tuned to reproduce the response of this detector, to deter- 
mine the best design for a larger-scale detector. 



The output of the Monte Carlo is the wavelength of 
photon hits on the PMT surfaces. In analysis, an energy- 
dependent quantum efficiency is applied, and a variable 
single photo-electron response is applied to each photo- 
electron. The trigger is modeled by requiring at least four 
PMT signals to be greater than a given threshold. The 
detector data is then scaled by the PE to ADC unit value 
determined from single photo-electron calibrations. 

There are three optical properties that can be tuned to 
reproduce the response in data: the wall reflectivity, the 
attenuation length of water, and, to a limited extent, the 
PMT quantum efficiency. The attenuation length of water 
in our detector was likely adversely affected by UV sta- 
bilizers in the polyethylene tank. The attenuation length 
from |9l [lOl [11] scaled down to an approximately 10 me- 
ter maximum gave results consistent with our data. The 
quantum efficiency is taken from Hamamatsu specifica- 
tions, and scaled down to account for losses in collection 
efficiency from stray magnetic fields: the total efficiency 
peaks at 22%. Finally, the wall reflectivity has a total 
reflectivity of 90% with a 5% specular component. 

Good agreement between data and Monte Carlo has 
been obtained with both AmBe source data inside the de- 
tector (Figure |9| an d the 252 Cf source data outside the 
detector (Figure 10). For the AmBe comparison, the cali- 



bration arm is included in the Monte Carlo. 



4. Discussion and Conclusions 

Based on the average number of neutrons produced 
from our 252 Cf source, as well as the average number of 
neutrons produced in reactor grade plutonium (RGP), we 
have determined that our current setup can detect approx- 
imately 150 g equivalent of RGP over background at two 
meter standoff from the detector edge in 20 seconds. The 
solid angle presented by our detector at this distance is 
2.4% of 4tt. We have assumed 8% fraction of 240 Pu in RGP 
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Figure 9: The pure neutron spectrum obtained from the tagged- 
calibration arm positioned inside the detector (black) compared with 
the MC spectrum from 5 MeV neutrons generated at the center of 
the detector (red). 
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Figure 10: The pure neutron spectrum obtained from a 252 Cf source 
positioned 7 inches from the detector edge (black) compared with 
the MC spectrum from 5 MeV neutrons, 7 inches from the detector 
edge. 

(87,000 n/s/kg). Figure [TT| is the signal over background 
for each cut in Table [I] as a function of radial distance 
from the detector. At one meter from the detector edge, 
or 6.4% of 47r, we have a signal to background ratio of 
0.547 ± 0.120 after all cuts have been made. Before cuts, 
the signal over background was 0.101 ± 0.00810. At two 
meters, the signal over background before cuts is 0.00560 
± 0.00755 and increases to 0.113 ± 0.0927 after all cuts 
have been made. Since our detector is intended as a cor- 
related event detector, it is expected that high solid angle 
coverage is required for efficient operation. However, with 
the proper event-level cuts our source is detectable even 
at two meters away: without requiring neutron-neutron 
pairs, the signal over background is 0.0702 ± 0.0154. 
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Figure 11: The signal-over-background ratio for 252 Cf adjacent to 
the detector, one meter from the detector's edge, and two meters 
from the detector edge. The colors correspond to no cuts (black), 
the charge cut on the current event (blue), the current charge and 
muon veto cut (purple), the charge cut on the current and previ- 
ous event (orange) resulting in neutron-neutron pairs, the neutron- 
neutron pairs plus the timing cut (red), and the neutron- neutron 
pairs plus the timing cut with the muon veto (green). There is an 
arbitrary offset in the x-axis as a visual aid: the green curve corre- 
sponds to the true radial position. 

4.1. Next Generation Detector 

There are many ways to improve the performance of 
this detection technology. Aside from optical and geomet- 
rical optimizations, decreasing the neutron capture win- 
dow by increasing the concentration of GdCl3 will result 
in a decrease in the rate of accidental coincidences. An ex- 
amination of this detector's 252 Cf calibration data shows 
that if the neutron capture window was decreased from 35 
/is to 10 /is, and assuming the same number of neutrons are 
captured, the percentage of accidental coincidences drops 
from 53% of all events surviving the current and last event 
energy cuts to 23%. In other words, if all neutrons are cap- 
tured within 10 /is, then for every 4.3 coincident events, 
one would be accidental. This is a significant improvement 
compared to our current concentration in which all neu- 
trons are captured within 35/is: for every 1.9 coincident 
events, one is accidental. It has been shown that there is 
no measurable impact on water attenuation length at the 
level of 0.2% GdCi3 [4]. An increase in GdCi3 concentra- 
tion to the 0.3% level is needed to obtain a characteristic 
capture time of ~10 /is. Further study may be warranted 
in this area. 

4.2. Conclusions 

We have successfully operated a large-scale gadolinium- 
doped water Cerenkov detector and characterized its neu- 
tron detection performance for the purposes of SNM mon- 
itoring. Using a tagged americium beryllium source, the 
raw neutron detection efficiency has been measured in the 
center of the detector at 70% and outside the detector 
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at 31%. Event-level cuts have been established to max- 
imize the detection of correlated pairs of neutrons emit- 
ted simultaneously from a 252 Cf fission source. We have 
demonstrated detection of approximately 150 g equivalent 
of RGP over background within 20 seconds for a solid angle 
coverage of 2.4% of 47r; an optimized geometry is expected 
to perform even better. 

The neutron capture response of the detector has been 
reproduced in Geant4 for two different source positions, 
using the water attenuation length, the PMT quantum 
efficiency, and the wall reflectivity as tuning parameters. 
The tuned value of each parameter is within the expected 
range for the water quality, PMTs, and Tyvek reflectivity. 
Future work will use this model to design an optimized 
radiation portal monitoring system for detection of corre- 
lated neutrons from undeclared fission sources within cargo 
containers. 
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